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The  Air  Force  Weapons  (Laboratory  ( AFWL )4 performed  Hydrodynamics  Unlimited 
(HULL)  calculations  ofl  the  air  blast  over  a  dam  for  two  yields  and  two 
pressure  regions.  A  £4-fth  calculation  included  a  rigid  blockhouse  at  the 
foot  of  the  dam.  Although  the  shielding  effect  of  the  dam  reduced  the 
overpressureyfrf  thp  incident,  ♦■hy  reflection  of  the  shock  from 

the  valley  floor  raised  the  peak  overpressure  up  to  at  least  40%  of  the  free 
air  value.  In  almost  every  case,  the  overpressure  impulses  near  the  foot 
of  the  dam  were  greater  than  or  equal  to  the  free  air  values.  The  rigid  ^ 
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SECTION  I 
INTRODUCTION 


a0at"„:TeS  °7he°reUCa’  ca,“,ati™s  “f  "“Clear  weapon  air  blasts  prop- 

g  ver  a  dam  structure  has  been  completed.  The  bursts  were  located  over 
water  behind  the  dam  at  such  distances  to  give  free  field  „»i 

figure  i,  1  t  *  yneS/$q  ™  <5°  a"d  10  PS1')  leyond  the  dam. 

F  re  1  illustrates  the  problem  configuration.  These  calculations  were  done 

With  weapon  yields  of  50  and  1000  KT  (209.15  and  4183.0  terajoules). 

the  cal  2'D  CarteSian  hydr0d»na,"1<:s  code  (ref.  1)  was  used  to  perform 

th  calculations.  For  computational  convenience  and  cost-effectiveness,  the 

calculations  were  simplified  by  specifying  a  problem  mesh  to  include  only  the 

around  the  dam  and  using  a  left  boundary  condition  that  replicates  a  blast 

wave  entering  from  the  left.  one  of  the  calculations  a  perfectly  "g"d 

structure  was  inserted  to  Investigate  the  loading  of  the  blast  wave. 
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SECTION  II 

COMPUTATIONAL  TECHNIQUES 

Two  computer  codes  were  combined  in  these  theoretical  calculations.  The 
Air  Force  Weapons  Laboratory  (AFWL)  HULL  (Hydrodynamics  Unlimited)  and  the  LAMB 
(Low  Altitude  Multiburst)  codes  were  used.  Both  codes  were  developed  at  the 
AFWL  during  the  last  5  years.  The  HULL  code  is  a  hydrodynamics  code  that  solves 
the  hyperbolic  equations  by  a  two-step  finite  difference  technique.  On  the  other 
hand,  LAMB  is  a  system  code  that  models  nuclear  phenomenology. 

1 .  HULL  CODE 

The  HULL  two-dimensional  version  was  used  to  predict  the  hydrodynamic  motion 
of  the  air  blast  waves  in  these  problems.  Basically  the  code  solves  the  follow¬ 
ing  equations. 

Conservation  of  Mass 

^  •  ft  =  0  (1) 

Conservations  of  Momentum 

p  +  =  "p9  ^ 

Conservation  of  Energy 


p  +  $  •  (PtF)  =  -pti  •  g 


(3) 


Equation  of  State 


P  =  P(p,E) 


(4) 


where 

P  =  pressure  in  dynes/cm2 
E  =  total  energy  in  ergs/gm 
0  =  fluid  velocity  in  cm/ sec 
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g  =  acceleration  of  gravity  in  cm/sec2 
p  =  material  density  in  g/cm3 
t,  T  =  time  in  seconds 

The  first  three  equations  are  approximated  by  a  finite  difference  technique 
developed  by  Matuska*.  There  are  two  phases  in  the  solution;  phase  1  solves 
equations  (2)  and  (3)  and  phase  2  solves  equation  (1)  by  fluxing  mass  across 

cell  boundaries. 

For  the  case  of  a  two-dimensional  problem,  Cartesian  coordinates  may  be 
used.  The  first  three  equations  are  then  rewritten  as 


where 

R  =  horizontal  coordinate 
z  =  vertical  coordinate 

u  =  component  of  velocity  in  radial  direction 
V  =  component  of  velocity  in  vertical  direction 

To  establish  the  finite  difference  analogs  to  equations  (5)  through  (8),  we 
consider  a  discrete  subset  of  F(R,Z,T)  by  defining 

F ( I » J  ,N)  =  F(R(1)  ,Z(J) *T(N ) ) 

where  R(I),  Z(J),  and  T (N )  are  particular  values  of  R,Z,  and  T,  respectively, 
and  the  I,  J,  and  N  assume  integer  values  in  the  range  1  to  I MAX  for  I,  1  to 
JMAX  for  J,  and  0  to  UMAX  for  N.  The  R(I)  and  Z(J)  are  defined  in  terms  of  a 

given  set  of  DR(I)  and  DZ(J)  such  that 
*Matuska,  Dan,  Private  Communication,  1975. 
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R( I )  =  R(0)+(SUM,K=1 ,  I- 1 , (DR(K) ) )+DR( I)/2  FOR  1=2 , . . .  , IMAX 

R(l)=R(0)+DR(l)/2 

Z(J )=Z(0)+(SUM,K=1,J-1 ,(DZ(K) ) )+DZ(J)/2  FOR  J-2.....JHAX 

Z( l)=Z(0)+DZ(l)/2 

where  R(O)  and  Z(O)  have  some  specified  values. 

The  hydrodynamic  variables  RHO  (material  density,  p),  U,  V,  and  E  are 
defined  for  each  set  of  coordinates  (I,J)  at  a  particular  time  T(N).  The  pres¬ 
sure  P(I,J,N)  is  defined  at  each  point  by  the  equation  of  state  (equation  4). 

Interpolated  values  of  the  hydrodynamic  variables  of  the  form  F(I+1/2,J,N), 
F(I,J+1/2,N),  or  F(I,J,N+l/2),  or  similar  forms,  are  defined  in  terms  of  the 
F(I ,J,N).  In  general 

F(I+1/2,J  ,N)  =  (F(I  +  l,J,i  ,+F(I  ,J,N))/2 
and 


F(I,J+1/2,N)  =  ( F( I  ,J+1 ,N)+F(I,J,N))/2. 

This  defintion  will  apply  except  where  explicitly  noted. 

a.  Phase  I 

Using  the  above  convention,  we  can  write  the  finite  difference  analogs  to 
equations  (6)  through  (8)  as 

U(I, J,N+1)«U(I,J,N)-DT*(P( 1  +  1/2,  J,N+l/2)-P( 1-1/2, J.N+1/2)) 

/(RH0(I,J,N)*DR(I))  (9) 

V(I,J,N+l)=V(I,J,N)-DT*(P(I,J+l/2,N+l/2)-P(I,J-l/2,N+l/2)) 

/(RHO (I ,J ,N )*DZ (J ) )-DT*G( J)  (10) 

E(I,J,N+l)=E(I,J,N)-DT/RH0(I,J,N)*((P(I+l/2, J.N+1/2) 

*U( 1+1/2, J,N+l/2)-P(I-l/2, J.N+1/2) 

*U( 1-1/2 , J ,N+l/2))/DR( I)+(P(I ,J+l/2,N+l/2) 

*V( I ,J+l/2,N+l/2)-P( I,J-1/2,N+1/2)*V( I , J-1/2.N+1/2) ) 
/DZ(J)-DT*V( I ,J,N+1)*G(J)  (11) 
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where 

DR(I)  =  horizontal  dimension  of  the  Ith  column. 

DZ(J)  =  vertical  dimension  of  the  Jth  row. 

DT  =  T(N+1 )-  T(N) 

R( 1  +  1/2 )  =  R(I)+DR(I)/2 

R( I -1/2 )  =  R ( I ) - 1/2  DR ( I ) 

Z( J+l/2 )  =  Z(J)+DZ(J)/2 

Z ( J- 1/2  )  =  Z ( J ) - 1/ 2DZ ( J ) 

G(J)  =  value  of  the  gravitational  constant  at  Z(J). 

All  the  values  appearing  in  equations  (9)  through  (11)  are  immediately 
known  except  the  time  advanced  (N+l/2)  values  for  pressure  and  velocity.  These 
time  advanced  values  are  used  so  that  the  approximations  to  the  partial  deriv¬ 
atives  appearing  in  equations  (6)  through  (8)  may  be  centered  in  time  and  space. 
In  the  case  where 

DR(I)  =  constant  for  I  =  1,...,IMAX 
and 

DZ( J )  =  constant  for  J  =  1,...,JMAX, 

this  produces  a  fully  second  order  accurate  difference  method.  In  a  region  where 
the  DR ( I )  and  DZ(J)  are  not  constant  the  second  order  accuracy  is  lost.  This 
adversely  affects  the  stability  of  the  first  phase  calculation.  The  amount  of 
instability  which  may  be  obtained  is  related  to  the  magnitude  of  the  incremental 
changes  in  DR( I )  and  CZ(J). 

Most  of  the  computations  1/;  the  first  phase  are  expended  in  obtaining  the 
time  advanced  values  for  pressure  and  velocity.  The  time  advanced  velocities 
are  obtained  by  differencing  equations  (6)  and  (7)  as 

U(I  +  1/2,J, N+l/2)  =  U(  1+1/2 ,J , N ) - DT / ( 2*RH0( 1+1/2, J, N+l/2 ) ) 

*((P(I+1,J,N)-P(I,J,N))/(R(I+1)-R(I)))  (12) 

V(I,  J+l/2, N+l/2)  =  V(  I  ,J+l/2 ,N)-DT/(2*RH0( I , J+l/2, N+l/2 ) ) 

*(  (P( I ,J+1 ,N)-P( I ,J  ,N))/(Z( J+l)-Z( J) ) ) 

-G( J+l/2)*DT/2 
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where 

G( J+l/2)  =  (G(U)+G(J+1) )/2. 

The  time  advanced  densities  appearing  in  equations  (12)  and  (13)  are  obtained 
by  differencing  equation  (5)  as 

RH0( 1+1/2, J.N+1/2)  =  RH0( I+1/2.J ,N) *( 1-DT/ (2*R( 1+1/2 ))*( R( 1+1) 

*U(I+1,J,N)-R(I)*U(I,J,N))/(R(I+1)-R(I)))  (14) 

RHO ( I , J+l/2, N+l/2)  =  RHO ( I ,J+1/2,N)*( 1-DT /2*( V( I ,J+1,N)-V(I,J,N)) 

/(Z( J+l )-Z(J ) ) )  (15) 

where 


RH0( 1+1/2 ,J ,N)  =  (M(I,J,N)+N(I+1,J,N))/(PI*(R( 1+3/2 )**2 
-R( I - 1/2 ) **2 ) *DZ { J ) ) 

RH0(I, J+l/2, N)  =  (M(I,J,N)+N(I,J+l,N))/(PI*(R(I+l/2)**2 
-R( 1-1/2 )**2)*( DZ( J)+DZ(J+1) ) ) 

and  the  mass  associated  with  a  point  (I,J,N)  is  defined  by 

M(I,J,N)  =  RH0( I ,J ,N)*(PI*(R( I+l/2)**2-R( I-l/2)**2 )*DZ( J ) ) 

where 

PI  =  3.14159. ..and  R( 1+3/2)  =  R( ( I+l)+l/2) . 

The  time  advanced  pressure  appearing  in  equations  (9)  through  (11)  requires 
a  little  more  effort.  First  an  alternative  energy  equation  can  be  obtained  from 
equations  (2)  and  (3)  as 


+  P(V-ti)  =  0 


An  effective  y  can  be  defined  by 


Y  = 


1  + 


P_ 

Pi 


(16) 


(17) 
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We  will  assume  for  the  purposes  of  calculating  a  half  time  step  advanced  pres¬ 
sure,  which  in  turn  is  used  in  approximating  the  partial  derivatives  in 
equations  (9)  through  (11),  that  the  Lagrangian  derivative  with  respect  to  time 
of  gamma  is  small  and  can  be  ignored.  In  application,  it  is  only  required  that 
the  change  in  gamma  at  a  particular  point  be  small  over  a  time  of  DT/2. 

Taking  the  Lagrangian  derivative  with  respect  to  time  in  equation  (17)  and 
using  equations  (1)  and  (16),  we  can  write 


~  +  (GAMMA) (P)($.t)  =  0  (18) 

Equation  (18)  is  used  to  obtain  time  advanced  pressures  given  by 

P( 1+1/2, J.N+1/2)  =  P( 1+1/2 ,J ,N)*( 1-DT*GAMMA( I+1/2,J,N)*(R(I+1) 

*U(  1+1, J ,N)-R( I )*U( I ,J,N))/(2*R( 1+1/2 )*(R( 1+1) 

-R( I ) ) ) ) 

P ( I ,J+l/2,N+l/2)  =  P( I , J+l/2 ,N)*( 1-DT*GAMMA( I ,J+l/2 ,N) 

*(V( I ,J+1,N)-V( I ,J,N) )/(2*(Z(J+l)-Z(J) ) ) ) 

where  gamma  is  obtained  from  equation  (17)  as 

GAMMA( 1  +  1/2, J,N)  =  1+P( 1  +  1/2, J ,N)/ (RH0( 1+1/2, J  ,N)*I ( 1+1/2, J  ,N) ) 

GAMMA(I , J+l/2 ,N )  =  1+P( I ,  J+l/2 ,N)/(RH0( I , J+l/2, N)*I (I , J+1/2.N) ) . 

All  quantities  needed  to  solve  equations  (9),  (10),  and  (11)  are  now 
defined.  Solution  of  these  equations  will  complete  a  second  order  accurate 
Lagrangian  calculation.  The  next  step  would  normally  be  that  of  transporting 
mesh  vertices.  Instead  we  choose  to  flux  the  hydrodynamic  quantities  to  retain 
the  orginal  mesh  configuration.  This  calculation  is  done  in  Phase  II. 

b.  Phase  II 

Changes  In  density  are  computed  in  phase  II  by  calculating  a  mass  flux  be¬ 
tween  mesh  points  and  then  transporting  the  appropriate  amount  of  mass  from  point 
to  point.  The  transported  mass  carries  with  it  a  proportionate  amount  of 
internal  energy  and  mt^ientum.  The  velocities  and  specific  internal  energy  are 
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then  redefined  at  each  mesh  point  by  conserving  momentum  and  total  energy  at 
that  point. 

The  mass  flux  between  mesh  points  is  defined  as  the  product  of  the  interpo¬ 
late  velocity,  the  density  as  defined  by  solution  of  equation  (i),  the  inter¬ 
mediate  cross  sectional  areai,  and  the  time  step. 

MF(I+1/2,J,N+1)  ■  Ut 1+1/2, J,N+l)*RH0(I+l/2,J,N+l)*2*PI*R(I+l/2) 

*I)Z(  J)*DT 

MF( I  ,J+l/2 ,N+1)  =  V(I, J+l/2, N+1)*RH0(I, J+l/2, N+1)*2*PI*R(I)*DR( I ) 

*DT 

where  the  time  advanced  densities  are  obtained  by  differencing  equation  (5)  as 

RH0(I+1/2,J,N+1)  =  RH0{ID,J,N)*(1-DT/R( I+1/2)*(R(I+1) 

*U(I+1,J,N+1)-R(I)*U(I,J,N+1))/(R(I+1)-R(I))) 

RH0( I ,  J+l/2, N+l )  «  RHO ( I , JD,N)*( 1-DT*(V( I , J+l ,N+1)-V( I , J,N+1) ) 
/(Z(J+1)-Z(J))) 

where 


ID  =  I 

IF 

U(I+1/2,J,N+1) 

GT 

=  1+1 

IF 

U(I+l/2, J.N+l) 

LT 

C_ 

o 

II 

c_ 

IF 

V( I , J+l/2 ,N+1) 

GT 

«  J+l 

IF 

V( I , J+l/2 ,N+1) 

LT 

This  is  the  classical  donor  cell  differencing  technique.  The  most  obvious 
advantages  of  this  technique  are  its  rigid  numerical  conservation  and  its 
stability.  This  scheme  also  insures  that  more  material  can  not  be  removed  from 
a  point  than  is  present. 

HULL  has  a  continuous  rezone  capability.  When  this  is  employed,  the 
interpolated  velocities  appearing  in  the  mass  flux  equations  are  replaced  by 

U( 1+1/2, J,N+1)-UR( 1+1/2, J,N+1) 

V(I,J+1/2,N+1)-VR(I ,J+1/2,N+1) 
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where  UR  and  VR  are  the  interpolated  grid  velocities  (determined  arbitrarily  by 
how  fast  one  wishes  to  transport  the  coordinate  grid). 

The  corresponding  momentum  fluxes  are 

UF(I+1/2,J,N+1)  =  MF(I+1/2,J,N+1)*U(ID,J,N+1) 

VF( 1+1/2, J, N+l)  =  MF( 1+1/2, J,N+1)*V( ID, J, N+l) 

UF( I , J+ 1/2 ,N+ 1 )  =  MF (I ,J+1/2,N+1)*U(I ,JD, N+ 1 ) 

VF ( I , J+ 1/2 , N+ 1 )  =  MF(I, J+l/2, N+1)*V(I,JD, N+l) 
and  the  energy  fluxes  are 

EF( 1+1/2, J, N+l)  *  MF( 1+1/2 ,J ,N+1)*E( ID, J ,N+1) 

EF( I , J , 1/2.N+1)  «  MF(I,J,1/2,N+1)*E(I,JD,N+1) 

When  these  quantities  are  fluxed,  final  values  for  mass,  density,  velocity, 
and  energy  are  computed  by 

M ( I , J )  «  M( I ,J ,N)+MF( 1-1/2 ,J ,N+1)+MF( I ,J-l/2 ,N+1) 

-MF( 1+1/2, J , N+l) -MF( I , J+l/2 , N+l) 

RH0(I,J)  =  M(I,J)/(PI*(R( I+l/2)**2-R( 1-1/2) *DZ(J ) ) 

U(I,J)  =  ( U ( I , J ,N+ 1 ) *M ( I ,J,N)+UF(I-1/2,J ,N+1)+UF( I , J-l/2 ,N+1) 

-UF( 1+1/2, J,N+1)-UF( I,  J+l/2,  N+l) )/M(I ,J) 

V ( I , J )  =  ( V( I ,J ,N+1 )*M( I ,J,N)+VF ( 1-1/2 ,J ,N+1)+VF( I , J- 1/2 ,N+1) 

-VF( 1 , J+l/2 ,N+1)-VF( 1+1/2, J ,N+1) )/M( I ,J) 

1(1, J)  =  (E(I,J,N+1)*M(I,J,N)+EF(I-1/2,J,N+1)+EF(I, J-l/2, N+l) 

-EF(I+1/2,J,N+1)-EF(I,J+1/2,N+1)-(U(I,J)**2+V(I,J)**2) 

*M(I,J)/2)/M(I,J) 

E  ( I ,  J )  =  I(I,J)+(U(I,J)**2+V(I,J)**2)/2 

where  the  lack  of  a  time  specification  indicates  final  values  for  this  time 
step. 
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2 .  LAMB  CODE 

The  LAMB  code  is  a  three-dimensional  system  code.  It  is  basically  a  model 
based  upon  two  and  three-dimensional  HULL  calculations  (ref.  2).  The  phe¬ 
nomenology  that  LAMB  can  predict  includes  shock-shock,  shock-fireball,  fireball 
fireball,  shock-ground  and  firebal 1 -ground  interactions.  Since  it  is  a  model, 
it  can  quickly  predict  a  complete  hydrodynamic  description  for  a  freefield 
environment  at  a  point  in  space  and  time. 


2.  Needham,  C.,  Matuska,  D.,  Bauer,  B.,  and  Whitaker,  W. ,  Air  Force  Weapons 
Laboratory  LAMB  Model ,  AFWL  Technical  Note,  1972. 
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SECTION  III 
PROCEDURE 

The  computational  method  that  has  been  employed  in  these  calculations  is  a 
hybrid  of  the  normal  hydrodynamic  calculations  with  the  LAMB  model.  Since  the 
LAMB  code  predicts  only  the  free  field  environment,  the  hydrodynamics  code  HULL 
must  -be  used  to  predict  the  fluid  flow  over  or  around  objects.  The  configuration 
of  the  problem  is  a  blast  wave  running  off  the  top  of  a  man-made  dam.  It  is  a 
two-dimensional  problem  that  can  be  solved  in  slab  geometry;  that  is,  two- 
dimensional  Cartesian  coordinates.  Initial  conditions  assume  a  standard  tem¬ 
perate  atmosphere  at  sea  level  with  respect  to  the  base  of  the  dam. 

1.  MESH  SIZE 

In  order  to  sharply  define  the  dam  and  the  structure  used  for  the  compu¬ 
tations,  the  size  of  the  rectangular  mesh  chosen  was  200  columns  and  130  rows. 

\  The  area  covered  by  the  mesh  is  indicated  in  figure  1.  In  addition,  the  area 

from  the  left  boundary  to  the  right  of  the  structure  was  finely  zoned  with  cells 
being  equal  to  1  meter  in  the  horizontal  and  vertical  directions.  Beyond  the 
structure,  the  cell  sizes  were  incrementally  increased  10  percent  in  the  hori¬ 
zontal  direction. 

2.  BOUNDARY  CONDITIONS 

There  are  two  types  of  boundary  conditions  normally  used  in  the  HULL  code. 

A  reflective  boundary  is  utilized  at  the  center  of  symmetry  or  whenever  a 
perfectly  reflecting  surface  is  used.  The  transmissive  boundary  utilized  is 
simply  a  set  of  mesh  points  constrained  to  be  ambient.  It  is  perhaps  better 
i  termed  an  observer  of  a  boundary. 

For  these  problems  we  have  used  reflective  boundaries  at  the  top  and  bottom 
and  transmissive  at  the  right.  At  the  left  boundary  the  time  dependent  hydro- 
dynamic  values  were  preset  from  LAMB  at  each  time  step.  The  AFWL  LAMB  code 
was  used  to  obtain  the  correct  variables  for  the  specified  weapon  yield  as  a 
function  of  time  for  the  boundary.  A  good  representation  for  the  air  blast 
i  waveform  can  be  obtained  in  this  manner. 

In  order  to  represent  the  dam  as  well  as  the  structure  within  the  problem 
mesh,  the  island  construct  of  the  HULL  code  was  used.  Islands  are  simply  a  name 
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for  an  algorithm  which  inserts  reflecting  boundaries  for  any  geometry  within  the 
mesh.  As  a  result,  the  dam  and  structure  below  the  dam  are  outlined  with  reflec¬ 
ting  boundaries  and  the  cells  within  these  boundaries  become  inactive. 

3.  STATIONS 

In  an  effort  to  monitor  the  loading  on  the  structure  and  the  physical  param¬ 
eters  of  the  shock  as  a  function  of  time,  stations  were  placed  at  various  points 
in  the  mesh.  T‘ e  hydrodynamic  variables  of  these  points  are  recorded  for  every 
time  step  and  stored  on  a  separate  data  tape.  Figure  2  indicates  the  location 
of  the  stations  in  the  mesh. 

From  the  station  data  one  can  compute  the  arrival  time  of  the  shock  front, 
the  peak  overpressure  and  its  impulse,  the  dynamic  pressure  and  the  dynamic  pres¬ 
sure  impulse. 
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SECTION  IV 
RESULTS 

This  study  of  air  blast  loading  included  four  combinations  of  yield  and 
range.  They  were  1  megaton  and  50  kilotons  (4183.0  and  209.15  terajoules)  at 
the  50  psi  and  10  psi  (3.45  x  10b  and  6.89  x  10s  dynes/sq  cm)  ranges.  In  'ddi- 
tion,  a  fifth  problem,  50.017,  included  a  perfectly  rigid  blockhousp  strut  re  at 
the  foot  of  the  front  face  of  the  dam.  Table  1  summarizes  the  calculations. 

Table  1 

SUMMARY  OF  PROBLEMS 


HULL 

PROBLEM 

YIELD 

BURST 

POSITION 

HEIGHT 

OF  BURST 

NOMINAL 

FREE  FIELD 
OVERPRESSURE 

NUMBER 

X0  (M) 

ABOVE  WATER 

dynes/sq  cm 

51.013 

1000 

1360.9 

0 

3.45  x  106  (50 

psi) 

51 .014 

1000 

3098.9  , 

0 

6.89  x  105  (10 

psi) 

51 .015 

50 

427.0 

0 

3.45  x  106  (50 

psi) 

51.016 

50 

1063.1 

0 

6.89  x  105  (70 

psi) 

51 .017* 

50 

427.0 

0 

3.45  x  106  (50 

psi) 

*This 

problem  includes  the  structure 

at  the  foot 

of  the  dam. 

To  understand  the  overall  phenomenology,  consider  the  two-dimensional  contour 
and  vector  plots  in  appendix  A  for  problems  51.015  and  51.017,  50  kilotons  at  the 
50  psi  range. 

For  the  first  30  meters  from  the  left  side  of  the  computational  mesh,  the 
blast  travels  over  the  top  of  the  dam  reservoir,  which  is  assumed  to  be  rigid. 
When  the  blast  reaches  the  front  of  the  dam,  it  expands  downward  reducing  the 
overpressure. 

The  reduction  in  peak  pressure  is  shown  by  the  contour  plots  of  appendix  A. 

At  a  time  of  480  milliseconds  for  problem  51.0170  the  peak  pressure  near  the  face 
of  the  dam  is  approximately  2.5  atmospheres  while  the  free  air  shock  is  over  3.4 
atmospheres.  As  the  shock  front  expands  down  the  face  of  the  dam,  a  rarefaction 
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begins  to  develop.  Examination  of  values  in  table  2  indicates  this  effect. 
Stations  24,  25,  and  26  are  at  the  same  ground  range  but  at  increasing  heights 
above  the  tcp  of  the  dam.  The  peak  overpressure  increases  with  the  height  above 
the  dam  from  2.77  x  10fi  dynes/sq  cm  at  a  point  even  with  the  dam  top  to  3.61  x 
10!  dynes/sq  cm  a  distance  32  m  above  the  dam.  The  impulses  are  similarly 
affected.  The  reduction  in  pressure  and  impulse  above  the  dam  is  a  result  of 
mass  and  energy  flow  down  the  face  of  the  dam. 

From  the  time  the  shock  enters  the  left  boundary  to  a  time  of  480  milli¬ 
seconds,  problems  51.017  and  51.015  (with  and  without  the  structure)  are  identi¬ 
cal.  As  the  shock  impinges  upon  the  structure,  however,  the  flow  is  considerably 
modified. 

Consider  first  the  less  complex  shock  reflections  as  the  blast  wave  traverses 
the  flat  surface  with  no  structure  (problem  51.015). 

In  problem  51.015,  the  blockhouse  was  composed  of  air.  Trace  particles 
located  in  the  region  of  the  blockhouse  moved  with  local  air  velocities  as  the 
blast  swept  by.  They  simulate  the  debris  from  a  very  fragile  structure. 

The  shock  reaches  the  bottom  of  the  dam  at  a  time  of  just  over  0.5  second; 
at  this  point  the  shock  is  much  weaker  (about  half  the  strength)  than  the  free 
air  shock  above  the  dam.  By  a  time  of  0.53  second  the  free  air  shock  has 
passed  out  of  the  region  of  the  blockhouse  and  the  shock  near  the  base  of  the 
dam  has  reflected.  The  peak  reflected  pressure  is  greater  than  the  free  air 
pressure,  but  it  must  be  remembered  that  the  vertical  flow  is  stagnated  and 
therefore  there  is  little  dynamic  pressure  associated  with  this  overpressure. 

An  interesting  side  effect  is  the  circulation  pattern  developing  along  the 
face  of  the  dam.  This  is  the  same  effect  that  would  be  seen  over  the  trailing 
edge  of  an  airplane  wing.  The  flow  over  the  dam  creates  a  partial  vacuum  along 
the  face  of  the  dam.  This  flow  and  reduced  pressure  have  little  effect  on  the 

peak  pressure  near  the  base  of  the  dam,  but  do  reduce  the  positive  duration  and 
impulse  near  the  base. 

As  the  reflecting  shock  travels  away  from  the  base  of  the  dam,  the  effect  of 
the  dam  in  reducing  peak  pressure  becomes  less  important.  This  is  seen  by  the 
increasing  peak  pressures  between  stations  1  and  27  (from  3.6  to  4.1  atmospheres). 
The  reflected  pressures  at  stations  1  and  27  are  somewhat  greater  than  in  the  free 
air  stream  as  described  above;  however,  as  the  distance  between  the  shock  and  the 
dam  base  increases,  the  shock  characteristics  approach  those  of  the  free  air 
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Table  2. 

CALCULATION  OF  1  MT  EFFECT  ON  DAM  STRUCTURE  AT  50  PSI  RANGE  PROB  51.013 
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regl0n  The  overall  effect  of  the  da*  becoams  negligible  after  the  chock  has 
traveled  approximately  four  dam  heights,  in  this  case  400m. 

The  above  descriptioh  ahd  effects  are  modified  when  the  structure  at  the  base 
of  the  dam  is  included.  Now  consider  the  two-dimensional  contour^ 
plots  for  problem  51.017  (appendix  A).  The  blockhouse  is  pe  io„'betwee„ 

hiast  is  more  complicated  than  in  the  previous  example.  In  9 

te  dal  face  and  the  blockhouse,  the  shock  is  intensified  by  -  convergence, 

-  two  surfaces.  The  shock  strikes  the  upper  con,  of  the  strut-  J* 

imately  4*  ms.  The  reflected  pressu  e  at  th .  ■  ^  ^  ^  free  afr 

planar  surface  reflection,  as  expected, 

Shock  pressure.  -  The  center  of  the  roof  and "  cls^erably 
struck  at  very  nearly  the  same  time.  9  atmospheres).  The 

greater  than  the  initial  shoe  ™  ^  ls  completely  stagnated, 

shock  continues  down  the  front  face  o  h  structure,  and 

horizontally  and  vertically,  by  the  constraints  of  th e  d am,  t  e  t 

the  ground.  This  results  in  a  peak  up  the 

incident  wave  (9.0  vs.  1.6  atmospheres)  The  "  m1dpoint.  ^ 

rr:  “r,rrr...„ ... .  -  » 

percent  greater  than  the  initial  shock,  and  arrives  approximately  50  ms 
the  first. 

When  the  reflected  shock  reaches  ^  than 

pressure  has  dropped  to  about  150  percen  o 
the  free  air  incident  pressure. 

On  the  back  of  the  structure,  pressures  are  significantly  reduced,  and  at 
the  midpoint  the  peak  overpressure  is  less  than  1  atmosphere. 

..  ctrurture  tends  to  translate  the  structure. 

The  pressure  gradient  across  the  structure  tenos 

this  case  the  pressure  difference  is  about  4.3  atmospheres. 

Tables  2  and  3  list  the  peak  values  of  the  station  data.  Problem  51 .01  is 

■  .ho  in  nsi  (3  45  X  10‘  dynes/sq.cm.)  region.  From  the 
nominally  1  megaton  in  the  50  psi  (3.45  x  y  one  ca„ 

overpressure  at  station  26,  the  station  closest  to  the  free 

see  that  the  overpressure  is  7  percent  below  the  nominal  value. 

The  station  data  show  that  the  effect  of  the  dam  is  to  th, i  incident 

shock  by  rapid  expansion.  The  Incident  overpressure  is  one  half  of 
va°ue  (at  station  26).  However,  the  reflection  of  the  shock  off  the  ground 
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Indicates  that  positive  phase  is  incomplete 
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below  the  dam  increases  the  overpressures  nn 

value.  erpressures  on  the  ground  above  the  free  air 

Overpressure  in  the  reflected  a^ 

the  shock  expands  upward.  Overpressure  lucre"”  '"th  lnCreaS’n9  aU'tudf''  « 

the  reinforcement  of  the  ref  e  e  „d  ?  1nCr*aSl"9  ™9e'  f™ 

not  from  the  Mach  stem  wh  c  ,nC,deBt  Enforcement  Is 

the  blockhouse.  The  reflects  IT™  T  ^  the  "“nity 

strengthens  the  fl  C,°"  *  da™'  ^ 

portion  of  the  flow  away  from  the  dam. 

Although  the  peak  overpressure  is  usually  lower  rear  the  hi  u 
the  free  air  above  the  dam  the  *  th  blockhouse  than  in 

blockhouse  than  above  the  dam.  VGrPreSSUre  impulSe  is  near  the 

ohiectTLcTas';!:;:  ;:;vTssure  and 

dynamic  pressure.  The  latter  l's  reduced"  «eT^"  *  **  dra9  frC"" 

the  dam.  The  dynamic  pressure  impulses  were  also  red  d  re9’°n  fr°nt  °f 
only  enhancement  of  dynamic  pressures  occurred  t"e rert^^reT"' 

656  3re  Sma"  C°"lpared  t0  the  fred  a''n  horizontal  dynamic  pressure. ’9"’  “* 

Each  of  the  five  problems  differed  sliahtiw  f 
50  and  ,0  ps  1  (3.45  x  ,0‘  and  6.89  x  ,o>  dynes/s,  cm) ^  ° f 

ences  in  the  inputs,  the  calculations  showed  g  a  eelT  Th7o 
phenomenology  was  the  same  fnr  thQ  k,  9  e  ent-  The  overall 

the  dam  protected  the  blockhouse  more  VlL  "cat  eTof  However’ 

of  50  psi.  ®  PS1  than  for  the  cases 
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APPENDIX  A 
SNAPSHOTS  IN  TIME 

The  following  information  will  be  useful  in  interpreting  the  three  types  of 
plots  that  follow. 

1.  Vector  Velocity  Plots:  These  plots  show  only  the  largest  velocities.  The 
scale  is  in  kilometers/second  and  the  magnitude  is  based  on  the  scale  in  the 
upper  right  hand  corner  of  the  plots.  The  point  of  each  measurement  is  the 
beginning  (tail)  of  the  arrow. 

2.  Relative  Pressure  Contour  Plots:  Relative  pressure  is  the  ratio  of  the  local 
pressure  to  the  ambient  atmospheric  pressure.  The  contour  lines  connect  points 
of  equal  relative  pressure  and  are  analogous  to  isobars.  The  contour  values 
given  in  the  upper  right  hand  comer  are  dimensionless.  D  x  1  is  the  width  of 
the  first  computational  column,  and  min  and  max  refer  to  the  minimum  and  maximum 
relative  pressure  values.  Note  that  the  dots  in  the  various  plots  represent 
different  things  in  the  same  plot. 

a.  Some  of  the  dots  merely  denote  measurement  stations  as  can  be  seen  by 
referring  to  figure  2.  These  dots  do  not  move  with  the  calculation. 

b.  Rigid  structures  are  outlined  by  dots  for  convenience.  The  dam  and 
structure  are  thus  both  outlined  by  dots.  These  structures  do  not  move  as  the 
calculation  progresses. 

3.  Trace  particles  have  been  placed  where  the  structure  would  have  been  in 
some  of  the  calculations.  These  trace  particles  move  exactly  with  the  flow  and 
"trace"  the  movement  of  any  fluid  displacement. 
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APPENDIX  B 
STATION  PLOTS 

Station  plots  are  the  calculational  equivalents  of  experimental  sensors.  At 
the  prescribed  station  points,  all  available  calculated  data  are  recorded.  The 
plots  show  the  value  of  the  hydrodynamic  variables  with  respect  to  time.  Note 
that  the  station  number  with  its  location  can  be  determined  by  referring  to 
figure  2  or  by  the  coordinates  given  in  the  upper  right  hand  corner  of  the  plot. 
The  radius  is  the  X  coordinate  and  the  height  is  the  Y  coordinate  with  respect 
to  the  lower  left  hand  corner  of  the  computational  mesh. 
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DtSTKIBI  TION  (Continued) 


CO 
UK  L 

ATTN.  AMXBR-TB 
ATTN:  AMXBH-BL 
ATTN:  AMXBR-RL 


CO 

Plcatinny  Arsenal 

ATTN:  SMUPA-TN 


CO 

liSAHO 
Ch.  Eng 

Dept,  of  the  Army 

ATTN:  DAKN-RDM 


Dir. 

USA  Eng.  WWES 
ATTN:  WESUL 
ATTN:  WESSS,  Joo  Zelasko 


Dir, 

NHL 

ATTN:  2027 

Dir. 

NRL/EOTPO 
ATTN:  5503 

CO/Dir. 

NEL 

ATTN:  4223 

CDR 

NWC 

ATTN:  753 

osc 

Dept,  of  the  Navy 
CO 

NWEF 

ATTN:  ADS 
DDR&E 

ATTN:  Asst.  Dir.  ,  Strat.  Wpns. 

Dir. 

DIA 

ATTN:  DI-7D 
ATTN:  DI-3 

Dir.  OSD,  ARPA 
ATTN:  NMR 

CDR 
FC  DNA 

ATTN:  FCPR 

ESD 

ATTN:  DEE 

Dir. 

WSEGp 

ATTN:  Doc.  Con. 

JSTPS 

ATTN:  JLTW 


DISTRIBUTION  (Continued) 
ERDA 

ATTN:  Lib. 


LLL 

ATTN:  Lib. 


DDC 

2  cy  ATTN:  TC 

GE  TEMPO 

ATTN:  T.  Barret 

NHL 

ATTN:  Code  2070 
ATTN:  J.  Boris 

NSWC 

ATTN:  Code  730 
ATTN:  J.  Peters 

ATTN:  Exp.  Eff.  Div.  ,  Mr.  Hornig 

DNA 

2  cy  ATTN:  DDST,  Dr.  Atkins,  Mr.  P.  Haas 

3  cy  ATTN'  STTL 

6  cy  ATTN:  SPSS,  E.  Sevln,  Captain  J.  Stockton, 

Major  T.  Strong,  K.  Goering  &  Dr.  Ulrich 

4  cy  ATTN:  SPAS,  J.  Moulton,  M.  Rubenstein, 

Captain  Garrison  &  Major  Anderson 
ATTN:  STAP 

3  cy  ATTN:  RAAE,  P.  Fleming,  Captain  Mueller  & 
Major  Bfgcni 

ATTN:  RATN,  CDR  Aldorson 
ARPA 

ATTN:  Lioutenant  Colonel  Whitaker 
JSTPS 

ATTN :  Maj  G  recne 

DIA 

ATTN:  Mr.  Castlcbury 

Sandla  Lab. 

ATTN:  Org.  3141 

ATTN:  Div.  1111 

ATTN:  Org.  1723,  C.  E.  Barr 

ATTN:  Div.  5644,  J.  Reed 

3  cy  ATTN:  M.  Merrit,  L.  Vortman  &  D.  Dahlgrcn 

Sandla  Lab, 

ATTN:  Org.  8000,  Dr.  Cook 
ATTN:  R.  Scholer 
ATTN:  T,  Gold 
ATTN:  J.  Mansfield 

LLL 

ATTN:  L-51,  M.  Hanson 

ATTN:  TID 

ATTN:  Jeff  Thompson 

LASL 

ATTN:  Rpt,  Lib, 

2  cy  ATTN:  J-10,  E.  Jonos  &  R.  Gentry 

Aerospace  Corp, 

ATTN:  Lib. 

ATTN:  R.  L.  Beok 
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DISTHIBUTION  <C’ontinuod) 

Bnj  Eng.  Co. 

ATTN:  V.  Kenner,  MS-44 
ATTN:  C.  Wayne 
ATTN:  J.  Dobkins 

Me  Don,  Doug. 

8  cy  ATTN:  K.  Stone,  P  Lewis,  D,  Bonn,  |\  Thomas, 
D.  Hildebrand,  J.  Logan,  R.  ",aeh.  6 
H.  Hordman 

Kaman  Sci. 

ATTN:  D.  Pirio 
ATTN:  D.  Sachs 

Kaman  Av. 

ATTN:  J.  Reutenik 

MartMar 

4  ey  ATTN:  R.  Heffner,  C.  Washburn,  II  Suguichl  fe 
D.  Cooke 

MRC 

2  ey  ATTN;  R.  Christian  6  D.  Sowie 
SAI 

2  ey  ATTN:  C.  Busch  &  R.  Lowen 

SAI 

ATTN:  U.  Hove 

SAI 

ATTN:  H.  Hitiendahl 
SAI 

ATTN:  M.  Totrinen 
SAI 

ATTN:  H.  Ciemens 

BTL 

ATTfJ:  W.  Troutman 
Caispan 

ATTN:  R.  Delibcris 
ISI 

ATTN:  W.  Dudzink 
Lockheed 
MITRE 

2  cy  ATTN:  S.  Morin  6  J.  Brown 
CRT 

ATTN:  M.  Rosenblatt 
H  &  D  Assoc. 

4  cy  ATTN:  H.  Brode,  C.  Knowles,  J.  Carpenter  & 

R.  Lelevler 

MARA 

ATTN:  S.  Kahalas 

TRW  Sys,  Gp, 

ATTN:  Ub. 

2  cy  ATTN:  A.  Kuhl 

AT'i  N :  J.  Peterson 

AVCO 

ATTN:  G.  Grant 


DISTRIBUTION  (Continued) 

IDA 

ASD 

ATTN:  16  L 

ARM  for  Stf,  Col 
ATTN:  Lib. 

Ind.  Col.  Arm.  For. 

Nat.  War  Col. 

ATTN:  Lib. 

USAMC 

ATTN:  AMCHD-HP-B 
USAEC 

ATTN;  AMSEL-RD 
USNCEL 

ATTN:  Code  L31 
ATTN:  CECOfcr. 

MIT 

ATTN:  Dr.  Hansen 
Princeton  U. 

ATTN:  PPL,  Dr.  Bleakncy 
Un.  MI 

ATTN:  Rsch.  Soe.  Ofc. 

ATTN:  IS6T,  G.  Frantti 

St.  Louis  Un. 

ATTN:  Dr.  Kissiinger 

Boeing 

ATTN:  H.  Carlson 
GATC 

ATTN:  MRD  Div.  Ub. 

II6NSp.  Pro.  Div. 

ATTN:  S.  Smith 

URS 

SWR1 

ATTN:  M  Whitfieid  “ 

Un.  IL. 

ATTN:  Dr.  Newmark,  Tal.  Iatb. 
IITRI 
RAC 

ATTN:  Lib. 

ADC 

ATTN:  ADSWO 
SACLO 
T  AC  LOS 
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niSTKim  TIONJC  ontinued) 
a  rw  i , 

ATTN:  IK),  l)r.  Miner 


2  ev  ATTN 

Sl!L 

ATTN 

DE 

ATI  N 

DEV 

ATTN 

l)EX 

ATTN 

DY 

ATTN 

DYS 

ATTN 

DYV 

10  cy  ATTN 

DYT 

ATTN 

DYX 

20  cy  ATTN 

DYT, 

Capt  M. 

Stueker 

20  cy  ATTN 

1)YT, 

Cnpt  M. 

Fry 

ATTN 

SA 

ATTN 

SA1) 

ATTN 

SAS 

2  cy ATTN 

SAT, 

M  r  Sha 

rp  *  Mr.  Murphy 

ATTN:  SAW 

ATTN:  Charles  E.  Needham 
ATTN:  Burton  S.  Chambers,  III 
ATTN:  Maj  Gary  P.  Ganonfi 


Official  Hecord  Copy,  Capt  Mark  Fry/DYT 


